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ABSTSACT 


The interstellar diff-usion of galactic- civilizations is reexamined 
by potential theory; both numerical and analytical solutions are 
derived for' the non-linear partial differential equations which 
specify a range of relevant models, drawn from blast wave physics, 
soil science, and, especially, -population biology. An essential 
feature of these models' iS' that, for all civilizations, population 
growth must be limited by the carrying capacity of the environment. 
Dispersal is fundamentally a diffusion process; a density- dependent 
diffusivity describes' -interstellar emigration. We concentrate on- 
two models, the first describing zero population growth (ZPG) and 
the second which also includes local growth and saturation of 
a planetary population, and for which we find an asjdnptotic travel- . 
ling wave solution-. For both models the colonization wavefront’ 
expands slowly and uniformly but only the frontier worlds^ are 
sources of - further expansion. For nonlinear diffusion with growth 
and saturation, the colonization wavefront- from the nearest independ- 
ently-arisen galactic civilization can have reached the Earth only 
if its lifetime exceeds’ 3 x 10^ yrs . For -ZPG the corresponding 
number Is 1.3 x 10 yrs. ' For these numerical results depend on 
our choices for the specific emigration rate,, the distribution of 
colonrzable worlds, and, in the second model, the population growth 
rate; but the dependences on these- parameters ia entrancingly weak. 
We conclude that the- Earth is uncolonized not because interstellar 
spacefaring societies are rare, but because .there are too many 
worlds to be conquered in the plausible lifetimes of the colonization 



phase of nearby galactic civilizations. However, the radar and 
television announcement of an emerging technical civilization on 
Earth may cause a rapid response by nearby civilizations, presently 
motivated to reach our system directly, rather than by random 
diffusion. 



Alexander wept when he heard from 
Anaxarchns chat there was an infinite 
number of worlds; and his friends 
asking him if any accident had befallen 
him, he returned this answer: "Do you 

not think, it a matter worthy .of lamenta- 
tion that when there is such a vast 
multitude of them, we have not yet con- 
quered one?" 

• Plutarch , On the Tranquility 
o'f the Mind 

"Where are they?" is a famous and possibly apo- 
cryphal rquestion posed by Enrico Fermi at Los Alamos in 
the late 1940s . In a galaxy with 'v<10^^ stars with planets 
apparently abundant and the origin of life seemingly re- ■ 
quiring very general cosmic circumstances; with the selec- 
tive advantage of intelligence and technology obvious and 
with billions of years available for evolution, should not 
extraterrestrial intailigence be readily detectable? This, 
question is most often phrased in its as.trophysical context 
(see, for example, Shklovskii and Sagan, 1966; Dyson, 1966): 
should not extraterrestrial civilizations millions or, perhaps, 
even billions of years mors advanced then our own be capable 
of so altering cosmic objects and energy sources, that a range 
of astrophysicai phenomena should exist which cannot be un- 
derstood apart -from the. hypothesis -of intelligent origin? 

ORIGSiNSAL PAGE SS 
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An answer is somecimes offered in the negative (Shklovskii, 

1976); i.e.j the apparent absence of astrophysical phenomena 
of intelligent origin is taken to demonstrate that no 
very advanced civilizations exist — either because ther-e 
are as yet xmdetermined impediments to the evolution of 
technical civilizations or because such civilizations in- 
evitably self-destruct early in their histories. Alternatively^ 
it can be argued (Sagan, 1974) that there are a wide range of 
poorly understood astrophysical phenomena, some involving ex- 
tremely high energies, and that (Sagan, 1973a) the manifesta- 
tions of very advanced civilizations would be no more apparent 
to us than the design and function of human engineering arti- 
facts are to ants crawling upon their surfaces. Any astro- 
engineering activity that is so 

wasteful of energy as to be observable with our limited technolo 
might, by its very nature, be necessarily shortlived. The fact 
that we have not yet acquired compelling evidence of such 
activity since our rise as a technological society -- . although 
there are certainly unexplained major astrophysical phenomena 
is not evidence for the absence of extraterrestrial civiliza- 
tions. It may rather be that we were not looking in the right 
place, with the right instruments at the right time. In any 
case, the sole productive approach is the open-minded pursuit of 
physical explanations of astrophysical phenomena, where only a 
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serious failure might be evidence for extraterrestrial in- 
telligence, but where, in our present ignorance, the possi- 
bility remains entirely moot. 

However, there is another and more modest side to 

the question. Pioneers 10 and 11 and Voyagers 1 and 2 show 

that even our infant technical civilization is capable of in- 

-4 -1 

terstallar spaceflight, although at velocities '^^10 pc yr 
Accordingly, .should not civilizations only a little more ad- 
vanced than we be effortlessly plying the spaces between the 
stars CSagan, 1963 j Shklovskii and Sagan, 1966)? And, if so, 
should we not on Earth today have some evidence of interstellar 
visits? Occasional serious attempts to deal with this question 
always conclude that no persuasive evidence of past visits exists 
in human legends or artifacts (Shklovskii and Sagan, 1966; Sagan, 
1979); and the most widely touted claims of such visits uni- 
formly have another and more plausible explanation -- ' generally 
in the area of archaeology or hoax (see, for example, Story, 
1976).. Likewise, the extraterrestrial hypothesis of UFOs 
that is, the contention that we are today being visited — 
also exhibits, despite a great deal of study, no persuasive 
evidence (Sagan and Page, 1972; Klass, 1974). These discussions 
by no means exclude past or present visits to the Earth; they 
merely stress the absence of strong evidence for such visits. 
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I£ the presence of extraterrestrial civilizations were 
to imply visits to Earth then the absence of snch visits would 
imply the absence of such civilizations. This argument was 
put forward by Hart (1975) who advanced it primarily In the 
context of interstellar colonization, not exploration. He 
argues that even if there are energetic difficulties or safety 
hazards (for example, the induced cosmic ray flxix) to relati- 
vistic interstellar spaceflight, there should be no serious 
obstacles to interstellar flight at wO.l c, in which 
case the Galaxy would be traversed in <10 yr. Other 
possible objections to large scale interstellar spacefaring 
such as loss of motivation or self-destruction -- Hart 
argues are unlikely because they must apply to all galactic 
civilizations to explain the absence of extraterrestrials on 
Earth . 


We believe Hart's analysis is flawed on a number of 
counts. On the one hand there are universal social impedi- 
ments to cosmic imperialism of a sort which should apply to 
every galactic civilization. Von Hoemer (1973) has remarked 
that even interstellar colonization at the speed of light 
cannot solve the presenr human population explosion on the 
planet Earth, With our present exponential growth rate of 
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Y = 0,02 yr “ and the colonisation sphere expanding at the 
speed of Light, in 500 years the expansion voltme will 
have a radius of 50 pc with all habitable planets in that 
volume reaching the Earth’s present population density; 
thereafter the population growth rate must decline. This 
is, of course, an extreme example but it demonstrates that 
every society which is to avoid severe overcrowding and 
exhaustion of resources must practice 'stringent population 
control and actively maintain very small values of y. Since, 
in the long run, exponentials defeat power laws on -every planet 
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in the Galaxy -- independent of local biology, evolutionary 
history, and social customs — any analysis of this problem 
must consider low values of y and, in the limit, zero p.op.ula- . 
tion growth (ZPG) . A similar point has been made by Cox (1976) 
One of the objectives of the present paper is to explore in 
some detail the consequences of low values of y for inter- 
stellar colonization. 

Another possible social impediment to rapid inter- 
stellar colonization, with some conceivable claim to uni- 
versality, is immortality (Kevles, 1975). In extrapolating 
to advanced technological societies we often seem willing to 
imagine daring engineering developments such as interstellar 
spaceflight with velocities v = c, but make only the most 
modest extrapolations in biology or psychology. If aging is 
due, e.g., to the accumulation of somatic mutations, it is. 
conceivable that an advanced society will have essentially 
eliminated both disease and aging and that the limits to indivi 
dual longevity will then be set by accidents and by events in 
geological and stellar evolution. The motivations which we 
consider reasonable may seem very unnatural to such a society. 
To what extent, for example, is our motive for interstellar 
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colonization itself a quest for immortality -- which 
would be rendered unnecessary if personal immortality 
existed-? A population of immortal organisms, each leaving - 
a offspring every year, exhibits a growth rate identical 
to that of a population of organisms in which each lives 
for a single year and then dies leaving (a + offspring 
(Cole, 1954), But if there is a selective advantage for 
immortality, the relative number of offspring of immortals 
and mortals will adjust themselves in favor of the former. 

A society of immortals would be required to practice more 
stringent population control than a society of mortals. In 
addition, whatever its other charms, interstellar spaceflight 
must pose more serious hazards than residence on the home 
planet. To the extent that such predispositions are inherited, 
natural selection would tend in such a world to eliminate 
those individuals without a deep passion for the longest 
possible lifespans, assuming no differential replication be- 
tween those in favor of and those opposed to individual 
longe-vity. The net result might be^ a civilization with a 
profoimd commitment to stasis even on rather long cosmic time 
scales and a predisposition antithetical to' interstellar 
colonization. 
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With even modest annual growth rates in science 

6 8 

and technology it is clear that a civilization 10 or 10° 
years in advance of ours would have technological capabili- 
ties which for us would be indistinguishable from magirc. 

The possibility exists that such a civilization might impose 
an absolutely stringent galactic hegemony. The establishment 
of an unbreakable Codex Galactica imposing strict injunctions 
against colonization of or contact with already populated planets (Shklovskii’ 
and Sagan, 1966, p. 451) is by no means excluded. This idea, 
which has been called the Zoo Hypothesis by Bali (1973), is 
not irrelevant to our concerns merely because it- is by defini- 
tion unprovable; the detection,, for' example by' radio asmronomy, 
of civilizations on planets of other worlds would make the 
hypothesis ralsifiable. It is clearly impossible to prognosti- 
cate reliably on the social behavior of a hypothetical 
advanced civilization, but these examples are put forth to 
suggest -that there may be many, still less apparent, social 
impediments to extensive interstellar colonization. 

On the other hand. Hart seems to have made a fatal 
confusion ben^een the velocity of individual spacecraft and 
the velocity of a colonization wave front (see also Jones, 1976; Cox, 

197-6) , In the history of human explorarion and colonization 
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it is not a single exploratory society which launches all 
such ventures, but colonies of colonies and. higher order 
descendants. For example, it might be argued that American 
exploration of the solar system is a fifth-order Phoenician 
exploratory venture, the Phoenicians having settled Carthage, 
the Carthagenians having settled Iberian seaports, the 
Iberians having discovered America, and the United States 
exploring the solar system. Lebanon and Tunisia and Spain 
are not today spacafaring civilizations. There was a waiting 
time of centuries to milLenia before the colony acquired suffi- 
cient resources to initiate independent exploration or coloniza- 
tion, during which- period the parent ci v: :li?a rf n rr da eli -npd . If we 
imagine the nth order colonization of a. new planet, the time 
for acquiring an (n -r l)th order independent colonization 
capability will be substantially longer. 

Jones (1975, 1978) claims to show from a Monte Carlo 
calculation that even allowing for this colonization waiting 
time it is very difficult to make the colonization wave front move 
more slowly than about 10*!^ of the starship velocity — implying 
that a single expansionist power will colonize the galaxy in 

years if .the starship velocities are 0.1 c. Colonization • 
for Jonas provides an escape valve for an overpopulated civiliza- 
tion. Consequently, his colonization front progresses rapidly. 
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Since, in hia model, the number of colonists at any given 
time is proportional to the number of "saturated" planets, 
the asymptotic behavior of his model must be characterized 
by exponential population growth and by a steadily a_c.celera- 
ting colonization front. We believe this conclusion to be 
in error and below present population dynamics arguments on 
the velocity v of such interstellar colonial ventures. It 
is clear that a comprehensive approach to the problem must 
include considerations of the population growth rates; the 
related waiting time for the acquisition of an independent 
colonial capability on a newly colonized planet; the number 
of colonists which arrive on a virgin world with .each expedi- 
tion (the larger this number the smaller the waiting time 
until the next order colonial venture) ; the starship veloci- 
ties and search strategies; the abundance of untenanted 
planets available- for colonization; the lifetime of the 
colonial civilization and its descendants; and, not unrelated, 
the possibility of interaction between two independent ex- 
panding colonial wave fronts originating oh separate worlds. 

We wish to assess the assumption, implicit in pre- 
vious speculations on this subject, that the colonization 
wave front advances with an effective velocity comparable to 
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that of the individual starships. In our analysis, we 
adopt some of the methods conventionally employed in de- 
scribing growth processes and dispersal mechanisms in 
terrestrial ecosystems. Their mathematical development- 
provides a rigoroiis description of the stochastic 
properties of such processes and permits the use- of" poten- 
tial theory. Consequently, we can trace the evolution of' 
the distribution of the organis-ms that are of interest 
■without following the detailed behavior of each of its 
individtial members , 

The models we presently describe can be. used, to 
simulate the growth, saturation and expansion of an inter- 
stellar civilization ”in the mean.” Since the mathematical 
formulation of the problem is accomplished via potential 
theory, our models -are corpletely specified by a family of 
nonlinear partial differential eq-uations. The properties 
of these eq^uations can be examined analytically, as well as 
numerically, and the results are directly amenable to 
scaling. To make the discussion entirely clear it will be 
necessary -first-to review some elementary arguments. 
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In biological applications, these models have 
proved to be sufficiently general in describing all but 
the most catastrophic population dynamics phenomena. 
Moreover, the results they provide in the inter steldiar 
context are remarkably insensitive to the ‘detailed de- 
scription (within limits) of the model and the associated 
input parameters. Consequently, the solutions obtained 
below are, we believe, less vulnerable to charges of 
terrestrial chauvinism than many 'previous studies. 

We now review the population dynamics processes 
relevant to our discussion, explore various features -of 
their associated mathematical models, and discuss the im- 
plications of the model results. The most significant of 
these is that the expansion velocity of the colonization 
front is several orders of magnitude smaller than had been 
previously anticipated. Thus, the answer to the question 
"X'Jhere are they?” may well be that only now are they about 


to arrive. 
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1 . Hoiao'geneous Process es 

It is convenient, first, to consider the popnla- 
tion dynamics processes that occur in a homogeneous en- 
vironment. [In this instance, we ignore the factors that 
contribute to immigration and emigration, for example in 
humans. In animal species, biologists describe such move- 
ments as "dispersal," a term used to describe all factors 
that lead to the displacement -of an organism from its place 
of birth to the location at which it reproduces or dies . 
[See, for example, Howard (I960).] The simplest description 
of growth is that the time rate of change of a population 
V 'is proportional to the population, or 



( 1 ) 


where the growth rate, y, is the difference between the 
birth and death rates . Exponential population growth is 
often associated with Malthus ' name although the idea that 
populations increase geometrically seems to have evolved in 
the Middle Ages, if not in classical antiquity.—^ 


1/ Leonardo Pisano (Fibonacci) in the year 1202 at- 
tempted to reintroduce into Europe the study of algebra. He 
posed the problem, complicated by overlapping generations, of 
how rapidly rabbits would reproduce from an initial pair in 
the first month and then a second pair in the second month 
before becoming infertile. The ntimber of pairs capable of 
renroducing in a given month corresponds to the Fibonacci • 
■series 1, 2, 3; 5, 8, 13, , . . Asymptotically, the number 

of fertile rabbits increases by [(5^ + l)/2] each month. 
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Th.e human population groxv'th rate has remained 
relatively xmiform, with a numerical value of approximately 
0.056 yr from the da^;ra of civilization until the middle of 
the eighteenth century. Since that time, in a period 
representing only 0.02% of human history, the Earth has 
witnessed a tenfold increase in human numbers. Population 
experts have speculaced that, for the past two centuries, 
we have been undergoing a demographic transition . It 
began with a decline in the death rate precipitated by 
advances in public health and nutrition. Later, the birth 
rate declined, primarily because of changes in the 
perceived value of having children. In the interim, the 
growth rate mushroomed and only now shows signs of stabilizing. 
Malthus in his - Essay on Population , published in 1798, 
recognized the fundamental incomp atability between a geo- 
metrical progression in population and the finite capacity 
of the environment co support life. He argued that, in the absence 
of 'hDral restraint,” other forms of birth control then unknown or considered 
unconscionable, the combined ravages of war, famine and pestilenc 
would impose a limit to growth. 

A conceptually straightforward generalization of 


Cl) which is frequently employed in describing this situa- 
tion is 


( 2 ) 
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where f'(v) has Che properties 

f(0) = 1; df(y)/dv < 0; f(v^) = 0, (3) 

and where v describes the carrying capacity of the en- 
vironment as the maxirrum population t-hpt it csn 
support- The function f(v) decreases in a manner that 
reflects the diminishing ability of a particular environ- 
ment to -support population growth. A simple but representa- 
tive example of this function is 


f(v) = 


1 - v/v. 


V < V. 


V > Vs 


(4) 


.The corresponding differential equation for the population 
growth (the Pearl-Verhulst law) 




has a general solution given by 

v(t) = 


v-(0) Ug e"'''" 


[V, - v(0)] + v(0) 


(5) 


( 6 ) 


where v(0) is the initial value of v(t). This curve has a 
characteristic S shape and is referred to as the logistic. 
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This curve shares the qualitative properties comaon 
to virtually all models of homogeneous processes: exponen- 

tial (J^falthusian) growth for v << v , followed asymptoti- 
cally by saturation at 

We employ the Pearl-Verhuls t 
law extensively in what follows to describe homogeneous 
processes. A particularly readable account of the 
appropriate mathematical theory (as well as 
equations describing multispecies interaction or 
"predator-prey” systems) may be found in Davis (1962) . 

Coale (1974) provides an interesting introduction to the 
history of the human population. 

2. Heterogeneous Processes 

Many factors contribute to the rates of dispersal 
in animal species. Certain ones are innate; others are 

environmental in origin (Howard, 1960) . The simplest of 
these, random dispersal, is a mechanism that evolved in 
order to reduce the incidence of homozygotes . An animal 
such as the muskrat,' Ondatra zibethica L. , will head off in 
a randomly selected direction over some characteristic 

distance before it mates, reproduces and dies. Since in- 
breeding is in this way avoided, the frequency of genenic reconbination 
of deleterious recessive mutations is reduced. 
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F.rom a probabii-istic point of view, each muskrat 
executes (in the absence of population growth) a random 
walk or Brownian motion; On the other hand, the behavior 
of the muskrat- population distribution is described by. the 
diffusion equation. 

The equivalence- of Brownian motion and diffusion is 
well-known (see, for example, Chandrasekhar, 1943; or Karlin 
and Tay.lor, 1975). Since a firm imder standing of these, 
phenomena is essential for our discussion, and has been in- 
co-mpletely appreciated in some previous discussions of the 
problem, . we consider a simple and standard example -- the 
drunkard's walk. We have a linear array of lampposts, each a 
'distance 1 . from the next, that lights the street- on which both 
a bar and a drunkard's home are situated. The drunkard 
is able only to walk the distance between one. post and the 
next before he collapses. Ha then arises a time t later and 
sets' -out toward another lamppost. Ihe-itiira. interval bamveen consecutive 
drunkard- lamppost collisions is then t,- Unfortunately, . he has 
no sense of direction and it is equally likely that he will 
stagger to the -left as to the right. If is the proba- 
bility that he is. at lamppost n at time interval k, then the 
equation governing his behavior is readily observed to be 


k + 1 
P_ 


J * T 
n -r 1 



1 


( 7 ) 
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G'Hien posed in this manner, the drunkard’s walk is an 
example of a Markov process; eq. (7) is known as the 
Chapman-Kolmogorov relation.) Rewriting, we obtain' 


•P 


k + 1 
n 



P-. 


K 


n -I- 1 


P-n 


1 


T 2t r 

Taking the limit as / and x go to zero (but in a manner 
whereby x / 2 t remains fixed) and noting that, from a first 
and second order Taylor expansion, 


( 8 ) 


|£- (x = ni, t = kT) = 


k+1 k 
m ■ *^n " ^n 


3^p(x = t = kx) 


3x 


/im ^n + 1 


ok, k 
Zd- +0 -I 
• n - 1 


(9) 




(where we have explicitly made the transition from a dis- 
crete to a continuous system), we find eq, (8) becomes the 
diffusion equation 


ia 

3t 



where the diffusion coefficient D satisfies the relation 



(10) . 


( 11 ) 


Generally, X and x are called the collision mean free path 
and the collision time, respectively. The mathematical 
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limiting process just described cen be performed 
rigorously. The essence of this demonstration is that 
there e^sists a precise method for describing random, 
stochastic processes via- partial differential aquations, 
a method known as po.teutial theory. In practice, poten- 
tial theory affords a much simpler and more reliable means 
of determining the behavior of a given, population (whether it 
be composed of muskrats, drunkards or, we claim, extraterres- 
trials) than the associated random process (which, other- 
wise is only accessible to Monte Carlo methods) . 


Before returning to the specific problem- of 
matching the dispersal of muskrats , let us explore some 
properties of the diffusion .equation. As a generalisation 
of (IQ),:, the diffusion equation can be written 


3o_(x,t) ^ „ 

pt " L 


|D(x,t,p) 7p(5,t) 


( 12 ) 


where p(x,t) now describes the population density. Here, 
Che dimensionality of the process is increased, in order to 
best represent the problem at hand (cypically, terrestrial 
problems are two-dimensional) . Hera D, the diffusion' co- 
efficient, need not be constant but, by this formulation, 
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the conservation of the total population jD(x,t)dx is assured. 

(In such conservative cases, it is cotmon practice to employ 
the terms population density and probability density inter- 
changeably, although these quantities are of course' not equal. 

The ratio of the ‘first to the second is fixed, and “the pro- 
portionality constant is the total population.) By allo-wing 
the diffusion coefficient to vary, we can, for example, represent 
the influence posed by geographical or environmental factors. 

When muskrats face a mo'untain range, their advance is slowed. 

In the present model, it follows that they will travel over a 
shorter distance in a time x than would be the case for a flat 
plain; the diffusion coefficient drops correspondingly. As a 
result, the flow of emigration observed does not follow straight 
lines but describes an extremum path of least resistance. 

(Some sources tend to discount the utility of employing diffu- 
sive terms in population dynamics. They furnish the erroneous 
counterexample that Minnesota would have been colonized by 
Europeans long before California, if the diffusion equation 
were emnloyed. However, they overlook the point that D should 
be decreased to represent the difficulties faced in overland 
travel by settlers coming to Minnesota 

and increased - to describe the rela- 

tive ease of ocean travel to California.) We see later how a 
population- dependent diffusion coefficient can be employed to 
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dascribe social, and territorial influences on amigration. 

?or the sake of siarplicity in our later discussion of in- 
tarstailar colonization, we assume that space is homogeneous 
over a suitably large scale so that D will have no spatial, 
variation. We shall, howeven, permit the diffusion- coeffi- 
cient to be. population density dependent. 


Since we will later be concerned with the problem 

of assessing Che velocity of the colonization front, it is 

important to point, out that the velocity associated with 

the diffusion process described by eq,. (10) is not well- 

2 

dafined. For example, if we define <x > ,- using the 
probability density p (x, t) as 


<x^> 



P (x, t) 


dx 


(13) 


and integrate eq, (10) by parts- twice, we find 


^ <x^> = 2D (14) 

Thus, the q-uantity ^ <x'^>^, which might be considered 
representative of the velocity, varies as t. Einstein 
(1905), in his classic microscopic treatment of. Brownian 
motion, showed that, this anomaly only arose when timescales 
shorter than the collision time are. considered. He 
demonstrated rigorously that the correct velocity at small times 
is tl-. Another way of characterizing the velocity y is 
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to compare the continuity equation (where j = pv is the ■ 
flux) 

’lE+!-i=° C15) 

with (12); i.e., we make the association -(known as Pick's 
law) that 

j = -D Vp (16) 

or, alternatively, 

V = -D Vp (17) 

Thus, the velocity associated with a diffusion coefficient 
independent of population density varies inversely as the 
density e-folding length but has no dependence on the 
density itself. In. both of these descriptions, each 
measure of velocity possesses certain undesirable qualities 
which can be attributed to the fact that the solution to 
the diffusion equation (10) is spatially unconfined, because 
the Green's function 

G(x,t) “ (4'n'Bt)~‘^ exp Ex^/4Dt] (18) . 

causes any initial configuration to be propagated and 
redistributed instantaneously .over all space. We later 
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see how nonlinear diffusion can overcome this velocity 
■anomaly . 


In practical situations, the diffusion coefficient 
is determined using Fick's law, ,i. e., as the quotient of 
the flux and the density gradient. In theoretical con- 
texts , the diffusion coefficient must sometimes be deter- 
mined-using other considerations. For example, discretizing 
(10) -"^d using length and time scales Ax and At, respectively, 
we find 


-k + 1 


n 


= D — . 
Ax" 


^n + 1 


, k 

Pr, - 


n - 1 


k L • 2mt\ 

■ I?- 


(19) 


Thus,- the effective probability P of leaving site n for 
either site ,n - 1 or n + 1 in a. time At is given by 


P = 2DAt/Ax^ 


( 20 ) 


Alternatively, the diffusion coefficient in m dimensions 
may be determined using the probability P, viz. 

■D = PAx^/2mAt, (21) 


Despite oxrr misgivings about the linear diffusion 
equation, it Has proven remarkably successful in describing 
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the dispersal of muskrats, Skellam (1951), citing earlier 
work, shows how well this model depicts the spread of 
Ondatra z'lb'ethica L. since its introduction in central 
Europe in 1905. But muskrats are not extraterrestrials and 
the model must be adapted to include social influences on 
emigration patterns as well as the effects of population 
growth and saturation. (Skellam' s treatment corrects for 
population growth, -but only a posteriori ) . The possibility 
of an advanced civilization which practices zero population 
growth must be explored. Moreover, all this must be done 
in the context of estimating the velocity of an expanding 
population front. 

. 3 . Density- Pep en d ent Diffusion 
Consider, now, the case of a diffusion coefficient 
proportional to some power of the population density, 

D(p) = Dg (p/p^)^ 

For N>0, this generalization of the diffusion equation 
(12) has several conceptual advantages. From (17), we 
obtain a velocity that now has a population density de- 
pendence. From the probability expression (20) for the 
drunkard's walk, the likelihood that a drunkard will leave 


( 22 ) 
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for another place now depends on how many other drunkards 
are converging on his lamppost. This is precisely the 
feature we desire: the probability of mounting a coloniza- 

tion venture is a rapidly increasing function of density, 
even if emigration cannot relieve the population pressure. 

[A more pronounced density dependence would arise by using 
an exponential dependence on the density. Crank (1956) in- 
vestigated the problem numerically in certain isolated cases. 
The analytic properties of’ the solution, however, are not 
well-understood. For our purposes, the power law in (22) 
will suffice in describing a representative cross-section 
of density .dependences . ] 

For the moment, we continue to neglect 
the homogeneous processes of growth and saturation. This 
might not be an xmreasonable assumption in describing an 
advanced civilization. It seems likely that such a society 
would have undergone a long period of zero population growth 
while it developed the technological capability for inter- 
stellar flight. Otherwise, exponential growth would have 
diverted all available resources to the task of feeding' a 
burgeoning population. ZPG must be an essential ethic of 
any society which is able to colonize other planetary systems. 
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In the spherically syiranetric three dimensional 
case, the population-dependent diffusion equation becomes 


'iP 

3t 


'll 9 
_ — , 

9r 


D 


[_2_1 ^ ia 

3r 


oj 


( 23 ) 


In this fo3nn, the equation describes the ’’porous medi;mi 
problem” and is frequently encoimtered in hydrology and soil 
science. [In the latter application, one considers the 
flow of a fluid through a porous medium. The mass of the 
fluid is conserved and it is assumed to obey a poly tropic 
equation of state. Instead of the Euler force equation, 
the fluid is said to obey Darcy’s law which provides for a 
velocity that is proportional to the pressure gradient. 

Muskat (1937) provides a definitive treatment of the problem, 
Philip (1970) offers an historical review of the subject 
including a treatment of selected topics in the theory of 
porous media, ] 


A formally identical equation is encountered in 
the physics of high- temperature gases, the "thermal wave" 
problem: 


9T _ 1 ■ 9 
9t r 9r 



a i — 
9r 
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(24) 
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Here, T is the temperature and the diffusion coefficient 
is temperature- dependent as -would .be the case in, for 
example, an ionized gas. Such a circumstance arises im- 
mediately after a high yield nuclear explosion , -when 
temperatures are sufficiently great that the velocity 
associated -with radiative energy transport far exceeds 
the somd velocity. Since diffusion is now much more 
efficient at higher temperatures, it quickly acts 
to, equilibrate the temperature distribution of the 
radiatively-heated gas. There is a thin boxindary 
layer separating the hot gas- from the -ambient atmosphere 
that has not yet been heated by the radiation. In this 
boxmdary layer, called the ’’thermal wave front," the 
temperature -undergoes a precipitous decline. Consequently, 
the thermal wave front is slow to conduct heat and the net 
temperature distribution changes slowly. 

The overall result of these different effects is 
•that the temperature distribution preserves a constant 
shape. Qualitatively, the temperat-ure is relatively 
uniform, from the origin to the thermal wave front where 
it declines rapidly to zero . Quantita- 
tively, the temperature (expressed in units of the 
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temperature measured at the origin) is a function only 
of the radius (expressed as a fraction of the thermal 
wave front radius) o To illustrate this, consider the 
one- dimensional analogue of (24) , . .n 


H = _L J 

8t 9x \ 


aT 



From this equation, we see that the quantity 


CO 


Q = 



^oo 


T(x,t) dx 


is conserved (Q is proportional to the total thermal 
energy) . There is only one dimensionless combination 
of the coordinate x and the time t that can be obtained ■ 
in terms of a and 0 using (25) and (26) : 


^ ^ (aQ^t)l/ ^N+2) 

The quantity (Q^/at)^^ has the dimensions of 

temperature and a solution to (25) which preserves its 
shape is 


I 2\ l/(N + 2) 

T(x,t) = f( 5 ) 


The solution for f(?) (see Zel’dovich and Raizer, 1967, for 


details) is 


(25) 


(26) 


(27) 


( 28 ) 



>29- 


fC?) = 


1/Nr - 1/N- 

’O ^ ) 

J L o 


C. < 


(29) 


where 


= 


L + M 1 + r - -jN/(N + 2) 

+ ^ • r(%- --H T/T^) 

- r (1/N) (30) 


For the case N = 0 [l.e., eq. (10) taking D 5 a] , 


f(^) = (4 tt) ^ e 


-% ,-1^4 


The normalization eii^Jloyed provides 


fCO d? = 1 


(32) 


and, for N>0, the position of the thermal wave front, 
using (27) , is just 


=• (a qN + 


^f -^o 


(330 
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In Figure 1, fC?) is shown for N = 0, 1 and 2. From 
eqs. .C29) and C31)', we see that: 

a) if N = 0, the distribution is tinconfined; 

b) if 0<N<1, the distribution has a finite cut-off 
where the temperature gradient vanishes; 

c) if N = 1, the temperature distribution has a 
finite cut-off with a finite, non- vanishing • 
temperature gradient; and 

d) if N>1, the temperature distribution has a 
finite cut-off with an infinite temperature 
gradient. 

The particular curves given in Figure 1 correspond .to ,a 
Gaussian, a parabola and an ellipse, respectively. Apart 
from a differing normalization factor, the f(C) profiles 
for spherical symmetry [eq. (24)] are unchanged. This 
solution to the porous medium or thermal wave problem was 
obtained independently by 'Barenblatt (1952) and by Pattle 
(1959) , A thorough, physically motivated discussion of 
this problem is given in Zel’dovich and Raizer (1967). 

Rnerr (1977) reviews some of the mathematical questions that 
arise in the one- dimensional version of this problem. 
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As we have already noted, the shape of the 
teinperat'ure distribution is preserved in this solution. 

Such a solution is called "self-similar [In the 
Russian literature, solutions of this type are called 
"automodelled" solutions . Cotirant and Friedrichs 
(1943) were the first to consider such shape-preserving 
solutions in application to gas dynamic shocks ; 
they employed the term "progressive waves."] 

Generally speaking, similarity solutions are charac- 
terized by an equation of the type 

r, = At*^ (34) 

r 

where r- describes the radius of the phenomenon 
being studied (e.g., a shock front or a thermal wave), 
and A and a are constants , 

There are two types of self-similar solutions. In 
the first, a is obtained by dimensional analysis and A is 
selected in order to satisfy a physical conservation law. The 
thermal wave problem is a case in point and eqs,- (30) and 
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(33) display the relevant features of the scaling.-^ 

A second type of self- similar problem is charac- 
terized when a is selected to insure the existence of a 
mathematical solution. Since ct is not known from 
dimensional considerations, the- self-similar solution 
is usually obtained by exhaustive numerical tests. 

The standard example of such a problem is the implosion 
physics studied by Guderley (1942). In his case, a<0 and 
trial values of a were varied (while numerically integrating 
the fluid dynamics equations) until a value of a was found 
that satisfied his boundary conditions. A more pertinent 
example of Type II self- similar solutions is that asso- • 
ciated with the semilinear diffusion equation in one 
dimension, (This equation is often used 


Another example of such so-called Type I self- 
similar solutions is that obtained for a strong isotrotic 
blast wave, as when the thermal wave described earlier cools 
sufficiently for gas dynamic energy transport to be more 
efficient than diffusive dissipation. This solution was 
first obtained by Sedov (1959) in .an application to super- 
novae. He argued that the shocked gas ’would lose its memory 
of its "initial" conditions due to convection and would ap- 
proach a solution characterized only by the total energy S 
of the original blast and the mass density p of the unper- 
turbed gas. In that case, 

^f ~ ^f ^ 

describes the radius of the shock, where is a dimension- 
less number of order unity selected to assure that the blast 
energy equals the total internal and kinetic energy of the * 
shocked gas. It is not without some irony that the methods 
for treaning explosive processes, such as thermal' or blast 
waves, and, in particular, nucl^ear explosions, - are valuable 
in studying population dyhamics.'' 
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in population dynamics modelling'; we examine it in 
detail, shortly.) In that case, the 'solution for the 
density p(x,t) can be expres.sed as a "travelling wave" 
p (x - vt) where v is the velocity of the travelling wave 
C’uniform propagation regime" in the Russian literature). 
If we take x = x ^ X t =t , then the solution can 
be written as p(?) where 

Comparing this with (34) , we find that 



describes the position of the travelling wave front, 
where is a constant and v has the same role as a, 

I 

„ Although self -similar solutions are special 
solutions characteritad by a dependence on a single • 
dimensionless variable, 5, they appear to describe the 
"intermediate- asymptotic" behavior of a much broader class 
of initial, boundary and mixed problems. In practice (for 
example, in the case of the thermal wave), the solution 
often loses its memory of the initial and boundary condi- 
tions . Then, the solution evolves into one characterized 





(35) 


(36) 
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only by, e.g., conserved quantities Csuch. as 0) and 
parameters that specify the problem (such as a) . This 
property of seif-similar equations has been rigorously demon- 
strated only for the linear diffusion equation (10)-, -and the 

semilinear diffusion equation. Numerical and experimental 
results (where available) display the convergence of solu- 
tions to self-similar form under a broad array of initial 
conditions. Thus, self- similarity is not of interest 
simply because it facilitates the calculations, but also 
because it describes a seemingly genuine asymptotic feature 
of the solution (see Barenblatt and Zel4ovich, 1972) . 


IS 


To recapitulate, we observe that a population- 
dependent diffusion coefficient overcomes many of the 
problems inherent in the conventional diffusion equation. 

From the self- similar solution, we observe that the distribu- 
tion is confined . Moreover, the velocity of advance of the . 
wavefront is no longer anomalous but can be calculated 
directly from eq. (33) , for one dimension, or (for N>0) from 




t) 


1 

3N+2" 


(37) 
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where 


«i = 



1 

_ / S 1 \ 


11 

3N+2 1 

3N + Z 

■ ■ rl? + h) 

i 

3n + 2 

2“ - 1 



li 



(38) 


in three dimensions CZel'dovich and Raizer, 1967). The 
Telocity monotbnically decreases as a function of time. 
CThis comes about since the central density decreases 
causing the diffusion rate to diminish accordingly.) 
Although this model ignores homogeneous processes, it is 
a useful description of situations where the overall 
population is conseorved (ZPG) and 

provides .useful insight into models that combine both 
homogeneous and heterogeneous processes. 


4. ■ Diffusion and Growth 

As a first step toward simulating the combined 
effects of honiogeneous and heterogeneous processes, con- 
sider the equation 


If - ap + a a 


3r 


( 39 ) 


This is similar- to an approach first proposed by Gurney 
and Nisbet (1975). By employing a population-dependent 
diffusion coefficient, they suggest that one can describe 
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0RtSir“3AL PAGE IS 
OF POOR QUALITY 


one non-random cause of dispersal, the "directed-motion" 
associated with territorial drives. (They confined their 
attention to the case N = 1, although they did allow for 
an additive third term describing random motion ; this 
they called the biased random motion model.) Amplifying 
on this theme, Gurtin and MacCarey (1977) proposed a method 
of solution for the one-dimensional analogue of (39) , 


3p _ ■ 3 ^ 1 a 

Tt- - PP + — y kp 


(40) 


After making the substitutions 


~ -ut 
p = e ^ p 


-1 


1 


T = [(a-l)u] (exp [(a-l)ut]-l 


(41) 


-h 


X* = x(ka)“^, N = a - 1, 


we obtain 





(42) 


We have already discussed the -solution to the latter 
equation [via eqs. (25) through (30)]. However, because 
of the scalings employed in (41‘) , we observe that 
p(x=0, t) must ultimately increase exponentially with 
time. To eliminate this anomaly, we must prevent the 
population from growing above the local carrying capacity. 
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.5 , ■ Hie SeiniTinear Diffusion Equation 
We' wish to describe the behavior of an organism 
with advantageous genes which undergoes random dispersal 
while increasing in population locally according to the 
Pearl -Verhulst logistic law.. In one dimension, this can 
be written 


f =YP (1-p/p,) 

where y is the local growth rate, p is the local carrying 

w 

capacity (or -population- density saturation level) of the 
environment and D is- the local diffusion coefficient 
(assumed here to be population- indep ehdent) . Further, it 
is assumed that y, p and D are constant and do not vary 
spatially or temporally. Making the transformations 


(43) 


p’ '= p/pg> . t’ =• yt 

(44) 

x’ = x(D/y)“'^ , 


eq , ( 43 ) b ec omes 


‘ie. 

3t 


P(l.- P) 



( 45 ) 


(where, for convenience, the primes have been dropped). 
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Fisher (1937) and Kolmogoroff , Petrovsky and 
Piscounoff (1937) independently proposed that eq. (45) 
has a right-going travelling wave solution. {The 
relationship of the travelling wave to self- s'imi 1 a r- ;• 
solutions was reviewed in the discussion leading to eqs, 
(35) and (36).] Assuming, then, that p depends only on 
x-vt in (43) , we firs.t transform the velocity v into 

V = v(Dy)'^ 

where v is a dimensionless'- quantity , and, employing (44) 
and (45) , obtain the ordinary differential equation 

-V = P (1 - p.) + P 

dx“ 

We wish now to find the range of v for which eq.- (47) 
possesses a solution. (Recall from oior earlier discussion 
that Type II similarity -s.olutions are characterized by a 
dimensionless constant, in this instance v, whose value is 
chosen in order to insure the existence of a mathematical 
solution.) Since we assTomed a right-going travelling wave, 
we should expect the density p to decrease monotonically 
in the direction of increasing x. Further, we assume that 

lim p(x) = 0 


(46) 


(47) 


lim 

X-+-CO 


p(x) = 1 


( 48 ) 
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i,e., that the density vanishes in the region iinaffected 
by the travelling wave, while the source of the organism 
Cinfinitely far off to th.e left) is supported at the 
carrying capacity. Since p is assumed monotonic in s, 
we can replace the derivatives in (47) according to 


dx ^ ’ 


q = - 


dp 

ds 


Cq is defined so that it is necessarily positive- valned) 
Then, (47) becomes 


vq = p (1 - p) + q 




Near the origin (i.e., p = 0), this equation has the 
asymptotic representation 

— 1 do 

vq = p, + q 3^ 

Since we expect that q vanishes when p = 0, consider 
a solution of the form 


q = b 


n 

P 


where b and n are real -valued positive constants . Sub- 
stituting into (51) , we find that n = 1 and 

vb = 1 + b^ 


C49) 


(50) 


(51) 


(52) 


(53) 
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Solving for b in terms of v, we obtain 

b = ^ (v^' - 4)^J /2 (54) 

Therefore, in order to guarantee that b be real-valued, 
the discriminant must be positive and 

V > 2 (55) 

CThe case v ^ -2 corresponds to a left-going travelling 
wave . ) 


Although Fisher and Kolmogoroff ‘et 'al. obtained 
this result, Kolmogoroff et went on to show that the. . 
solution to the partial differential equation (45), 
initialized as a step function, converges to a right-going 
travelling wave with the minimal velocity v = 2. Although 
a continuous spectrum of propagation velocities v A 2 is 
in principle possible, only the solution corresponding to 
the extreme point of the spectrum (v = 2) can be an asymp- 
totic solution to the partial differential equation (45) . 
Therefore, the velocity of the front of advance of the ad- 
vantageous geneSj in Fisher's formulation of the problem, 
'^converges" to 2(Dy)^. 
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Although this model from population genetics can- 
be made more general by employing a broader description 
iE situ growth processes than that provided by the 
Pearl-Verhulst law, : such models are characterized by 
random dispersal, i..e., a population- independent diffusion 
coefficient (see, for example, the review article by 
Hadeler, 1977)-. Gaaracteristically, they have travelling 
wave, solutions which, when normalized in a similar fashion, 
also give v =* 2. Thus, despite the presence of a linear 
diffusion term, the associated velocity of advance is well- 
dafined [as contrasted with the absence of a self-consistent 
velocity for the linear diffusion model in the absence of 
population growth, eq. (10)].' 

There is, however, one feature of this model 
(apart from- the population-independent diffusion coeffi- 
cient) that is undesirable. From eq. (54), we see that 
b = 1 when v is minimal and, from (49) and (52) , we find 

p (x) e e~^ 

when p is small,. Therefore , the population ' s travelling 
wave front is unconfined. [This, again parallels the result 


( 56 ). 
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obtained for the linear diffusion model (10). J 

This model has been employed in describing the 
interaction of several species. For example,’ in a 
homogeneous environment, the equations describing a' 
predator, denoted by p-j^, and its prey, denoted by P 2 , 
might be written 

.3 

TT 

9 p 2 

TF ^3P2 " ^4PlP2 

where c^, C 2 , c^ and c^ are positive constants. In the 
absence of prey, the number of predators will decay 
exponentially. On the other hand, in the absence of 
predators, the number of prey will grow exponentially. ■ 
When both predators and prey are present, the growth rate 
of predators increases while that of the prey decreases. 
This formulation of the predator-prey relationship is 
generally associated with the names of Lotka and Vol terra 
(see Davis, 1962, for a simple treatment of the mathe- 
matical theory) , and describes accurately a number of well 
observed ecosystems. The solutions to eqs. (57) describe 
closed orbits in the P;j_~P 2 plane and are akin to the solu- 
tion for a harmonic oscillator in that p-, and P 2 will 
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oscillate in a pattern strictly determined by initial 

conditions,. When eqs . (57) are augmented by linear diffu- 
2 2 2 2 

sion terms D 8 pj^/3x' and D 3 p2/8x , respectively, oirr 
■equations assume more of the character of the semilinear 
diffusion equation (43) . An example of the predator-prey 
relationship is seen in herbivoroxis copepods (zooplankton) 
which feed on phytoplankton. In performing' a linearized 
perturbation analysis for such marine planktonic communities, 
Steele (1974) noted that all finite wavelength perturba- 
tions' will be damped and the ecosystem will tend to spatial 
homogeneity. Nonlinearities, however, could conceivably 
destabilize the situation. Levin (1976) and' McMurtrie (1978) 
have reviewed the multi-species population dynamics problem. 

In the case just described, diffusion acts to 
stabilize the system (at least in a linearized representa- 
tion) . In other circumstances , diffusion can produce a 
travelling-wave similar, in certain respects, to that of 
the semilinear diffusion equation (43). This feat-ure was 
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■ 3/- 

employed by Noble (1974)— in devising a model for the 
Bubonic Plague of 1347. His equations are 


|i = K I S - yl + D I 
||-=-KIS + DV^S 

O t 


(58) 


where S and 1 are the densities of the susceptible and 
the infected populations, y is the disease’s mortality 
rate, and D is the diffusion coefficient. The coefficient 
K describes the rate at which the plague is transmitted 
locally. (There is, in fact, a third ’’species” implicitly 
included here, the density of individuals B -who have "con- 
tracted the disease and either perished or recovered. The 
’’mortality” rate y describes the rate at which this third 
population density grows, SB/9t = yB. The total population 


2/ Noble compared his model predictions with the 
information available on the Black Death (hanger, 1964) . 

His results appear to match what is known about the geographic 
advance of the disease, ‘ In commenting on the applicability 
of his model to other epidemiological problems, he suggested 
that the spread of certain mass sociopsychological phenomena 
(such as new religions) might be sufficiently similar to that 
of plagues to be worth investigating by similar methods. It 
is possible that the spread of UFO belief systems might be 
described in this way. 
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JC^ + S + B) d^x then remains constant,) Eq. (58) 

represents the rate of change of the number of infectives 

within a small area as the rate of transitions (KIS) from 

the susceptible population less the removal rate (due to 

mortality, --jil,. and dispersal D 9 I/9x ).- Meanwhile, the 

rate of change of the susceptible population within a 

small area is a net loss due to the transition to the in- 

2 2 

fected population (-KIS) and dispersal (D 9 S/9x-). Noble 
performed a numerical integration of the one- dimensional 
analogue of eq s . (58). His simulation quickly evolved into 
a pair of travelling waves, one describing the infectives 
and the other, the susceptibles . The qualitative shape of 
the two population distributions is depicted in 'Figure 2. 
•Thus, we see that semilinear diffusion acting in a multi- 
species problem can produce a set of travelling waves just 
as we observed in the single species case'. 

6. ■ Density -Dependent Diffusion with growth 
and' S'aturation 

We have already described the need for a population 
dynamic model that combines local growth and saturation 
processes with spatial dispersal mechanisms that are 
population-density dependent. This model, in one dimension. 
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may be expressed as 


3p 

3r 


= Y P (1 - ^ 


9x 


D ( 


p vN 9p 

0 ^' 3x 
s 


where the dxffiasron constant D xs scaled xn order to 

>- 

■ 4/ 

provide the diffusion coefficient when p = Ps*— From 
the transformations (44) , we obtain an equation analogous 
to (45) , namely 


(59) 


iP 

3t 


= p(i - p) + 4 


N 


i£ 

3x 


(60) 


We have investigated the properties of the solution to 
eq. (60) for yarious N by numerical means. The 
equation was integrated using a conventional four-point- 
explicit integration scheme (Richtmeyer and Morton, 1967) . 


4/ Density-dependent diffusion does not make emi- 
gration a response to population pressure alone; Since the 
fliix of the population is given by -D (p/pg)^^3p/3x, the 
dens'ity gradient plays an all-important role,_ Thus, for 
example, there is little net movement of a population that 
is approximately loniform in density, even if the density is 
nearly saturated. 
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The density distrihution, following Kolsogoroff ^ al . 

C1937) , was initialized, as a step function 

( 


pCs, t=0 



X < 0 

X > 0 




Since we ware also interested in the spherically sTumetric 
case, we mraerically integrated 


(61) 


If = Pd- - p) + 


3 

^ 


X p 



(62) 


using a stable differencing scheme developed by Eisen 
(1967) in applications to the spherically' symmetric diffu- 
sion equation. The initial conditions employed were 

1 • r < 

0 r > r^ 


p (r, t == 0) = •( 


(63) 


where-r^ is 'some arbitrary radiuS', As a test of the 
integration methods, the case of H = 0’ was included- In 
both geometries, for N = 0, 1 and 2, the numerical solu- 
tion evolved into a travelling wave, (The asymptotic 
correspondence between the one dimensional and spherically 
symmetric cases is nor surprising. The only difference, be- 
tween their, respective equations is (2/r)p^ 3p/3r. Assuming 
that the solution tends to the form p (r - vt) , this 
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tem makes a contribution of order t~^ , which at large 
tinHs is negligible. ) Since, for N = 0, the numerically 
evolved travelling- wave has a velocity very near 2, we 
consider the asymptotic convergence of the solutions (for 
N - 1 and 2) to travelling v?aves to be significant. Figure 3 
illustrates the density profile for these three cases. 

In order to facilitate a comparison between the wave fronts, 
rhe density profiles are constructed so that density is 
one-half at the origin. 

Several other interesting fearures of the solution 
emerge from the computer analysis. First, the velocity v 
of the travelling wave j^eq. (46)j decreases with in- 
creasing N, The values obtained for N - 1 and 2 were ap- 
proximately 0.70 and 0.45 respectively. (Higher values of 
N were also investiganed. Fluctuations in v made these 
cases difficult to estimate, but the decline of v with in- 
creasing N was unmistakable.) Secondly, the density profile 
near the point of vanishing density resembled the self- similar 
thermal x<rave profiles for N = 1 and 2 which appear in Fig. 2. 
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Consider, therefore, the problem of finding a 
travelling wave solution to (60) . Employing the methods 
used earlier in our discussion of the semilinear diffu- 
sion equation, we write the equation for p (x) as • 


-V 




cx 


- p(l - p) + 



Then, making the same assumptions here as were employed 
in obtaining (50) , we have 


(64) 


vq = p(l - p) + q ^ (pNq) 


(65) 


As before, 'we obtain an approximate solution as a power 
law- [eq. (52)], namely 


( 


q 


V p 


1 -N 


V + (v^ - 4) 


N > 0 
N = 0 


( 66 ) 


Near p = 1, the solution is given approximately by 




(v^ + 4)^ - V 


(1 - P) 


for all N. Integrating (67) , we see that the density 
exponentially decays to unity as we move inside the tra- 
velling wave. Near the front of the travelling wave, the 
density (for the N = 0 case) decays exponentially. 


( 67 ) 
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However, for N nonzero, the solution is of the form 
pCx) ^ in V (x^- - x)J 

where denotes a cut-off point to the 'right of which 
the density is identically zero. This behavior is essen- 
tially the same as we observed in the thermal wave problem 
for N>0 [see eq. (29) and subsequent discussion] and ex- 
plains the striking similarity observed in comparing the 
N = 1 and 2 curves in Figs . 1 and 2 (the thermal wave and 
travelling wave fronts, respectively). What has occurred 
is that, for N > 0, the growth and saturation teirm in eq. 
(64) has no influence on the shape of the density profile 
near the cut-off point arid a sblutiori' reminiscent of the 
self- similar thermal wave problem emerges. The growth and 
saturation terms, however, have a controlling influence on 
the determination of v. 


Upon finding this asymptotic behavior, we noted 
that (67) was an exact solution to (65) for N = 1 if we 
set V = 2'^, A straightforviard integration then gives 


P (x) = { 


1 - exp [2“'^(x - X )], x<x 


( 68 ) 


0 , 


X > X 


(69) 
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This profile corresponds , within the nnmerical linita- 
tions posed by the integration scheme, to the compnted 
solution and the predicted travelling wave velocity of 2 ^ cor- 
responds well with, the mmerical result of 0.70. Although we 
are not at present able to demonstrate the stability of this 
solution (in the sense of Kolmegoroff et al. ), the con- 
vergence observed in the nxmierical solution to the profile . 
given by eq. (69). is a relatively convincing demonstration 
of some form of asymptotic stability. However, unlike the case of the 
Fisher-Kolmo.goroff. problem, we have explicitly obtained a 

5 / 

closed form expression for the population density profile,— 


5/ Lax (1969), in presenting the ninth John von 
Neumann uecture of the Society for Industrial and Applied 
Mathematics remarked on how fascinated von Neiimann had be- 
come by the possibility that patterns disclosed by nnnaerical 
calculations, might reveal entirely unsuspected properties of 
solutions of nonlinear differential equations. A case in 
point _is_ the Kortweg-deVries equation which. describes long 
waves over water and some wave phenomena in plasma physics. 

In his lectiure. Lax went on to describe how fcruskal and 
Zabusky discovered the existence of certain solitary wave 
solutions after studying motion pictures of the computations. 
In 'our case, the computed solutions to eq, (59) pointed out 
the existence of an asymptotic travelling wave solution 
and, in one instance, led to a closed- form analytic expression 
for the profile and its associated travelling' wave velocity. 
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7 . S-ummary of Mathematical Approaches 

We have reviewed the advantages (and limitations) 
of potential theory in describing population dynamic 
processes. We have observed that dispersal is essentially 
a diffusion process whereas population growth (if it takes 
place) is restricted by the carrying capacity of the en^ 
vironment. A density- dependent diffusion coefficient was 
shown to be capable of describing the desirability of 
emigration (all other factors being equal). Essentially, 
two basic models emerged. One model described a condition 
of zero population growth. The properties of its solution 
may be obtained from the porous medium or thermal wave 
problems. The second model also included local growth and 
saturation. A numerical study of its properties disclosed 
the existence of an asymptotic travelling wave solution. 

Two features emerged from our population dynamic 
modelling. First, the velocity of the colonization wave 
front can be estimated by dimensional analysis using the 
physical and other parameters that influence the evolution 
of the civilization. The parameters that contribute to 
this problem include the growth rate y and the diffusion co- 
efficient D. The latter, in turn, is related by eq. (21) to 
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the effective probability P that an individual member of a 
civilization will embark on such a venture, the mean distance 
Ax between habitable solar systems and the effective time 
span At between successive colonization ventures (to be dis- 
cussed below) . Second, the dimensionless multiplying co- 
efficient that appears in the expression for the velocity 
(5 or V for the zero population growth and travelling wave 
solutions, respectively) is of order xmity for all models 
and the explicit dependence of the velocity on the physical 
parameters is characterized by power laws with fractional 
exponents. This manifests itself in velocities that are 
remarkably and encouragingly insensitive to the particular 
model and parameters employed. In the next section, we see 
how both models, applied to interstellar ' colonization, yield 
a much lower colonization wave velocity than had previously 
been predicted. 

Before continuing, we should point out why our 
models provide results that contrast so strongly with those 
of Jones (1976) . We recall that, in his model, all saturated 
solar systems were sources of colonists. In all of our 
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diffusion models that include population growth and 
saturation, the flux of colonists goes as 9p/3r 
(where we now employ spherical synnnetry) . From our 
computed solutions, the flux is significant only in_ a 
thin layer at the population wave front. In both Jones* 
and our models, colonists are the source of population 
growth (all growth has vanished from the saturated in- 
terior region). In our case, however, colonists are 
produced only near the surface of the sphere enclosing 
an interstellar civilization. We stress that this 
qualitative understanding of the difference between our 
models and that due to Jones evolved only after exploring 
the methods and assumptions implicit in current population 
dynamic models . It was only after performing an extensive 
analytic and mmerical investigation of the equations of 
population dynamics that we were able to distinguish the 
two most prominent features of population dynamic models: 
expansion due to directed diffusion from the colonization 
front and a velocity of expansion that is, at the same time, 
characterized by but relatively insensitive to the physical 
and sociological parametrers involved, 



In each case, the rate of population increase 
dv/dt is proportional to the number of colonists 9 present, i.e 

dv / dt - A6 (70) 

In our case, 9 is proportional to the 

surface area of the saturated sphere containing the ex- 
panding interstellar civilization. Since v varies- as the 
radius of the sphere, r^, to the third .power, o-ur model, pre- 
dicts 

9 - ; (71) 

Jones' model predicts 

9 = B’v . ' (72) 

Combining xeqs . -(70) and (71.) gives ■ • 

V <= t^ ■ . (73) 


Therefore, our model predicts that r^ « t and we have 
a ■ calonization sphere that is undergoing uniform expan- 
sion. On the other hand, Jones' model [eqs. (70) and (72)j 
yields 

V « exp (AB't) . . (74) 


Consequently, his modal gives an exponentially growing 
colonization front velocity. The critical difference 
between Jones' model and ours is that our model supplies 
colonists only from the periphery of the civilization. 
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Since we expect the core of the civilization to be 
saturated, almost all available resources would by definition 
be committed to the continued maintenance of the local eco- 
system and few, if any, resources would, be available for the 
construction of starships capable of transporting and sup- 
porting colonists across interstellar space. Interstellar 
spaceflight for a substantial community of colonists at 
v-vc-must for any civilization be enormously expensive in 
propulsion and shielding. At v<<c, the cost remains enormous 
because the long transit times require elaborate life-support 
systems. Contemporary colonies of vlO^ persons in cislunar 
space (O'Neill, 1975, 1978) are estimated optimistically to 
cost >>$10^^. Even with generous allowances for future tech- 
nological progress a larger colony designed to traverse parsecs 

safely might be argued to cost substantially more than the 

13 

'available resources of the planet Earth (v $10 gross planetary 
product) . It is not clear that for either velocity regime, 
colonization ventures to more distant targets would be cor- 
respondingly more expensive. Except for extremely advanced 
and long-lived technical societies, colonization is possible 
only to the nearest star systems. Thus colonization is 
initiated from the periphery of an empire for two reasons. 
First, the transit distances and times are shorter. Second, 



- 57 - 


materials are not in such short supply and the launching 
of colonization ventures would not place as -unbearable a 
strain on the outposts of the empire as on its central 
worlds . 


8 . Results 

We now collect some relevant results of the pre- 
ceding sections , apply n-umerical values, and address our- 
selves' to the central question of whether the Earth should 
have been visited during geological time by extraterrestrial 
civilizations . 
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Ttie steady state mmiber N of extant civilizations 
more advanced than our oxm in the Milky Way Galaxy can be 
written as 

■■ N = fL, 

where L. is the mean lifetime of galactic civilizations 
in years and f-is a factor which combines the rate of star 
formation, the fraction of stars with planetary systems, 
the number of ecologically suitable planets per such 
system, the probability of the origin of life on a given 
otherwise suitable planet, and the likelihood of the evolu- 
tion of intelligence and technical civilization (for ex- 
tensive discussions, see, for example, Shklovskii and Sagan, 
1966; Sagan, 1973b). A conventional estimate, which 
is, however, no more than a semi- informed guess, is f 
vlO'^yr However, in the case of colonization the 
probability of the emergence of intelligent life and 
technical civilization on a given world approach unity 
and this factor increases perhaps to f vlO yr 
and Sagan, 1966, P- ‘^•51). 


(75) 


(Shklovskii 
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Assxnning stars in th.e I'lilky Way Galaxy to have 
a mean separation of 1 parsec ■(== 3.26 light years), the mean 
distance between advanced civilizations is 

A - C2.5 X IO^Vn)^/^ =X2-.5 X 10^^/fL).^'^^ pc ■ . ‘ (76) 

We here asstme 250 billion stars in the Galaxy. Hence- 
forth, distances will be measured in parsecs and times 

in years . With a conventional value (op cit . ) for in- 

6 

dependently arisen civilizations, of W = 10 , A =* 25 pc - 82 l.y. 
Sq. (76) is based on a spherical distribution of stars in 
the Galaxy; but the. inaccuracies for distances larger’ than 
the thickness of the Galaxy (v IQO pc) is still small enough 
that we may neglect it in this . problem. 

For o\n: density -dependent diffusion problem with 
growth and saturation, we have found that: the velocity v 
of the colonization front can be represented by 

v=v (Dy)^''^ (77) 

where v is a dimensionless constant of order unity; D is 
a diffusion coefficient that, when miiltiplied by the 
population density gradient, yields the correct outirard 
fliox of population, as it nears saturation, and y is the 
population growth rats. For the human population on the 
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Earth today, y <^10 . In the centuries prior to 1750, 

y was 'v6 x 10~^ (Coale, 1974). In the late Pleistocene 

-7 -8 

y was much smaller still, perhaps 10" or. 10 

We have found jeq. (21)] from finite differences 
in the standard m-dimensional diffusion equation that the 
probability of diffusing to another planet Ax away in time 
At is P = 2mD At/Ax^. This is in fact the probability of 
the entire population emigrating to another solar system at 
a distance Ax over an effective time At. This is in stricc 
agreement with the behavior of a particle undergoing Brownian 
motion. However, since the diffusion equation describes the 
behavior of a distribution of particles, it is conceptually 
useful to think of P as the expecration value of the fraction 
of the population that will emigrate to another world at a 
distance Ax in a time At. In the sense of an expectation 
value, then, we define the specific emigration rate P by 

Y = P/At = 2mD/Ax^ . (78) 

This quantity describes the fraction of the population that 
emigrates per unit time to the next solar system. Here, At 
is the time interval between successive colonial ventures. 

If the transit time exceeds this time interval then we 
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substitute the transit time for At, Conversely, the problem 
is paced by the time interval between successive coloniza- 
tions for very fast transit times, so that even the extremely 
xmlikely contingency of faster-than-light travel (through^ 
multiply- connected space in the vicinity of black holes, say) 
does not alter this analysis significantly. 

j 

It is impossible to“make a very reliable estimate 
of ¥ for extraterrestrial civilizations.. Very high values 
should be stupifyingly expensive. Very low values cannot 
provide even short-term relief of population pressure at the 
periphery of the civilization. In the third century B.G. the 
Chin emperor approved a proposal of Hsu Fu to launch a 
colonization mission ’’with several thousand yoxmg men and 
maidens to go and l.ook for the abodes of the immortals hidden 
in the Eastern Ocean," (Needham,' 1971). The fleet of sailing ■ 
ships was never heard from again and the emperor complained 
bitterly about the cost. Over the following centuries a 
number of further expeditions to the Pacific Ocean to find 
an elixer of immortality were mustered, but none on the scale 
of the venture of Hsu Fu (o£_. cit . ) . All these expeditions’ 
may have perished; alternatively some- may have colonized 
Japan and have been the origin of the non-Ainu population 
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ORsGSi^AL PAGE IS 
OF. POOR QUALITY 

there. At this time the total population of China was 
'^^3 X 10^. The next taj or distant emigrations from East and 
Northeast Asia were the Mongol invasions and the explora- 
tion and trade by the Ming Navy in the Indian Ocean in the 
fifteenth century. These numbers give 'f vlO” yr”^. During 
the eighteenth century European colonization of North 
America, 'i' was v3 x 10"^yr~^ (Potter, 1965). Proposals 
have been made recently for the launching into Earth orbit 

of self-contained closed ecological systems with vio'^ in- 

_6 

habitants (O'Neill, 1975, 1978). This represents "'^lO of 
the population of the Earth. A very difficult undertaking in 
the next few centuries would be to launch one such Space 
City into interstellar space every century.- This gives 
^ 'vl0~®yr~^. We adopt 3 x lO^^yr"^ >¥> lO'^yr**^, with a 
bias to the smaller value for a younger interstellar 
civilization. 

The existence of abundant extra-solar planetary 
systems seems very plausible but is still undemonstrated. 
Presumably a civilization wishing to colonize the Earth 
comes from an Earth-like world; although^ since schemes have 
been proposed to terraform Venus and Mars — that is, to 
convert them from rather different environments to ones 
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which resemble nhe Earrh (Sagan, I960; Sagan, 1973.c) 
xeno forming the Earth into some rather different environ- 
ment might not be- beyond the capabilities of an interstellar 
spacefaring civilization. We note that such xenoforming' 
will, add to the waiting time for a colony to develop, an 
independent colonizing capability. Simple theoretical 

I 

models of the formation of planetary systems from solar 
nebulae suggest that a wide variety of planetary systems 
can be formed depending on, for example, the initial mass 
density distribution function of the nebula;, in some schemes 
terrestrial planets appear to be formed abundantly (Isaacman 
and Sagan, 1975) . In this case the distance between Earth- 
like planets is probably < 3 pc (see also Dole, 1964). Hart 
(1978,- 1979) has argued chat only an extremely unlikely, set 
of circumstances has preserved the Earth's environment 
particularly abundant liquid water -- from the opposing 
chreats of a runaway greenhouse effect and a global ice age. 
These arguments, ignore plausible variations in some climatic 
parameters and neglect negative feedback loops in terrestrial 
climatology; ws consider them unlikely. But if they are valid 
they increase Ax possibly to v300 pc,. Accordingly, we Xvill 
allow to range from 3 pc to 10 pc , with a preferred value 
or .10 pc , 
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If we adopt m = 3, Ax^ = 10 pc^, and ^ = 10“^yr”^ 
we find a typical value of the diffusion coefficient to be 
D = 2 X 10 ^ pc^ yr ^ — 5 X 10^^ cm^sec"^, an absolutely 
enormous diffusion coefficient by planetary atmospheres 
standards . 


We are now ready to examine the apparent absence 

of extraterrestrial colonies on the Earth. Let t_,. be the 

mm 

minimum time for the advanced technical civilization nearest 
us in space to reach us af a diffusion wavefront velocity v. 
Then v t_,. - A. If t ._ > L, there should be no such 

colonies, as observed. Consequently, there is a critical 
lifetime for a spacefaring interstellar civilization 


= t . = A/v 

c mrn 


(79) 


For the Earth to have been colonized, the colonial civiliza- 
tion must have a lifetime in excess of L^. 


Combining eqs . (76) -(79) we find 

- 3.3 X 10^ (f/0.1)“^/^[(Y/10"^) (f/10"®) (Ax^/10)] ^s. (80) 

For comparison, we combine eqs. (23), (24), (26), (76) and (78) 
to obtain the critical lifetime for a thermal wave front 
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with N - 2 and find 

= 1.3 X (31) 

This is the case of strictly observed zero population growth, 

2 

Our 'free parameters, f, y, Ax and ? have normalized standard 
values, already discussed, with distance measured in pc and 
time in years . 

We have now obbained, the critical 

lifetime L for the cases L <<y”^ (zero population growth) and 

L^>>y"’^ (nonline'ar population diffusion with growth and 

saturation). From eq, (80), we see that increases as y 

decreases. When y is set eqtial to L in (80), we obtain a 

9 ^ 

value for of 4.3 x 10 yr.s when f, •'? and Ax assume their 
normalized standard values. Although (80) is not strictly 
valid in this case, this value for suggests, as we in- 
tuitively expect, a smooth transition between our two models, 
i.e,, as -Y decreases, increases until it assumes the ZPG 
value . 


We see that for the colonization wavefront of the 
nearest technical civilization to have reached the Earth, 
the lifetime. of that civilization must exceed 30 million yrs . 
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for nonlinear diffusion with growth and saturation. Such 
a civilization will have been intensely occupied in the 
colonization of more than 200,000 planetary systems before 
reaching the Earth, some 64 pc away. Many colonial empires , 
of vast extent by terrestrial standards, may still occupy 
only an insignificant voliame of the Milky Way Galaxy, and 
not embrace the nearest independently arisen technical 
civilization whether it has colonial ambitions or not. 

If strict ZPG is observed, L is of the order of 

c 

the age of the ■ universe and it would be unlikely in the 
extreme for us to observe such a population. Variations in 
f would almost certainly increase the critical lifetime while 
variations in V will have an, almost negligible effect. 

Therefore, strict ZPG readily explains our failure to. observe 
nearby extraterrestrial civilizations and we will confine our 
remaining considerations to cases including population growth 
and saturation. 

Our conclusions for the case of density-dependent 
diffusion with population growth and saturation seem to be 
interestingly insensitive to the choice of input parameters. 
Reasonable variations in f, to allow for extensive colonization. 
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changes by a factor '>^3.. We believe thafy cannot be 
Giuch larger than lO^'^yr"^ to be consistent with the very 
powerful population pressures which any colonial empire 
must have experienced earlier in its history than the time 
of extensive interstellar colonisation, (In addition, one 
possible means of interstellar transport — the "genera- 
tion-ship” in which a much later generation of colonists 

arrives- than that which originally sat out — imposes an 

2 

extremely strict regime of Z?G on the settlers.) Ax might 

3 2 

be as large as 10 pc , meaning that the colonists must be 

in the practice of searching cotsparatively large volumes of space be 

fore establishing a colony. ¥ conceivably could be much 

larger than 10 yr" , although values as large as those 

typical for the- European colonization of Jlorth America in the 

eighteenth century seem prohibitive for interstellar space- 

2 3 2 14-1 

flight. But, even adopting Ax = 10 pc and 3 x 10 yr 

5 

we find that vlO vr. Even with massive interstellar 

C r 

colonization efforts, involving lO'X of the base population 
every century, the Earth would not have been visited unless 
the colonizing civilization, were very long-lived. 
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The effective velocity of the colonization wave- 
front is given by 

V = 1.2 X 10~^ [(y/IO”*^) (Y/10”®) (Ax^/10)] ^ pc ;yr"^ , 

Note from dimensional analysis that this is the only 

2 

combination of y, Ax , and. which will make a velocity, 

2 

With nominal values of y, Ax and H', the colonization wave- 

— fi — 1 

front has an effective velocity '^lO” pc yr~ 
considerably slower even than the interstellar velocities 
of the Pioneer 10 and 11 and Voyager 1 and 2 spacecraft. 


(82) 
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Givi liza-t ions with, liieciiaes ''^10 to LO yr, 
on the other hand, could have colonized the. Earth by now 
even i£ they originate some considerably greater distance 
from US' 'Within, the -Milky Way Galasy. However, for; reasons 
which we have already mentioned — dealing with the 
evolution of such very advanced societies -- we believe 
that their motivations for colonization may have altered ut- 
terly, and that their science and technology may be extremely 
different from anything •we- can recognize or even imagine. 

We have fottid that only technical civilizations with 
lifetimes' vio^ to 3 x 10^ years could, have initiated a 
colonial wavefront which has reached the Earth from the 
nearest technical civilization. The lifetimes of extra- 
terrestrial technical. civilizations are, of' course, highly 
uncertain but many v7orkers ‘in this field have hypothesized 
that only a tiny fraction of civilizations survive beyond 
10^ to 10^ years (see, for example, Shklovskii and Sagan, 
1966). In addition, there is the question of phase 
relations in the evolution of technical, civilizations. In 
the solar neighborhood there are_ unlikely to be many stars 
with habitable planets which are significantly older than 
the s-un, and we certainly expect' a substantial waiting 
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time -- perhaps several times 10 yr -- for the 
indigenous evolution of technical societies. There are 
unlikely to he any colonizing civilizations in the solar 
neighborhood >>3 x 10^ yr old. ■ 

It may be useful to divide technical civilizations 

6 6 

into two categories; young (<<10° yr -old) and old (>>10 
yr • old) . Old civilizations are unlikely to be colonial 
powers in the sense we are describing. Yoting civilizations 
may be embarked on extensive colonial ventures but, compared 
to the volume of the Milky Way, these are of very limited 
scope. There may be empires of tens or even hundreds of 
thousands of worlds. .But it is implausible that a Galactic 
Empire exists: there are simply too many worlds to conquer. 

The number Of habitable planets in the Galaxy may 
exceed the number of individuals on the Earth. Might the 
political evolution of the Galaxy to some extent parallel 
that of the Earth, in which individuals start out more or 
less alone, and eventually form larger family group-s., tribes, 
city/states, nations, super-powers, and, perhaps eventually, 
a single global state? In that case the Earth is at the 
present time still in the situation of a lone individual 
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one who even wonders if there are any' 'Other people at all. 

The stage we have been describing in this paper, when there 
are tens or himdreds of thousands ■ of colonial worlds for 
each interstellar civilization, corresponds to a much 
later stage in human history, roughly that of the- first 
city/ states. It is only at this point that groups interact. 

Th-us Star Wars, if there are any, occur at the level of 
city/ states — Athens versus Sparta. 

It is interesting to speculate on the likelihood 
of such. -encounters . For a random, uniform distribution 'of 
civilizations the probability p of’ an interaction would be- 
have as 

p < (d/A)^ (83) 

where:^A is the mean distance between civilizations ,|eq. (76^ 
and d is the -mean value of the maximum radius of an. interstellar 
empire, i.e,-, 

d = vL = v(Dy)^ I, . (84) 

Using (76), (79), (83) and (84) we find 

p < (L/L^)^ (85) 

for L < L . Thus, the likelihood of two civilizations 
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interacting is remote unless L 'v L^. (This argument de- 
pends of course on a given civilization growing in a 
manner that is independent of its neighbors.) 

If the probability of Star Wars goes as if 
the canonical lifetime for contact to be just occurring 

■7 

is 10 yrs; and if we were to imagine all civilizations at 
D A 

L = 10 yrs, only 1 in 10^ civilizations will be interacting. 

* 6 

Since N = 10°, “^lOO will be interacting at a given time. If 
5 8 

L = 10 yrs , then 1 in 10 civilizations will be interacting 

andj for N < 10 civilizations, there are no Star Wars. 

Thus, if warfare is an illness which is outgrown in the first 
5 • - 

10 years or so of the existence of a galactic civilization 
then there never are interstellar conflicts. 

■ The slow propagation speed of the colonization 
wavefront does not exclude specific exploratory missions 
to targets of particular astrophysical or biological 
interest. The wavefront of the nearest colonial civiliza- 
tion may be tens or even only a few pc from the Earth and 
able to visit our solar system in relatively short times in 
the likely case that starship velocities > 0.1 c exist. 

They have not visited the Earth .because there is no motivation, no apparait 
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way for a neighboring civilization to know of the existence 
of our developing interstellar civilization. This situation 
may change in the relatively near future. Military and 
astronomical radar systems and, especially, commercial tele- 
vision have been for the last several decades generating an 
immistakable radio signature, now vlO pc away, of the existence 
of an emergent technical society in the vicinity of the Sun. 
Depending on the disposition of nearby civilizations , their 
colonies and their exploratory vessels, and how carefully they 
are examining the radio spectrum in eavesdropping mode, the 
leakage radiation from our technology may stimulate an explora- 
tory mission from the nearest galactic civilization within the 
next century or less. 
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Figure 1. Thermal wave profiles for diffusion coeffi- 
cients with a power law density dependence. 

Figure 2. Wave of advance of the Bubonic Plague, 

Figure 3. Travelling wave front for three model equations. 




Figure 1 


Density 



Figure 2 




